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1. Introduction 



$_i ■ The notion of non-integer dimension although developed by mathematicians in 

late 19th - early 20th centuries has found its wide application in physics only in the 
last few decades. This application was mainly initiated by the pioneering work of B. 
Mandelbrot who attracted attention to the fact that "fractal" geometry provides an 
appropriate framework for the description of a whole range of complex structures 
that are formed from smaller subunits [ |I| . Whereas such structures are characterized 
by appropriate fractal dimensions, their growth and spatial correlations [ ||| are 
described by a (non-trivial) spectrum of multifractal (MF) dimensions [ |3j. [| . 

In this paper, we study the properties of diffusion phenomena in the presence of 
an absorbing polymer. This provides another example of a multifractal phenomenon 
in condensed matter physics. To this end we use the model proposed by Cates and 
Witten and derive the MF spectrum in the frames of a field theoretical formalism 
using the renormalization group [ |6], |7| method. The MF spectrum is related to the 
spectrum of exponents governing scaling properties of co-polymer stars [ 9. We 
calculate this spectrum to the third order of perturbation theory and report the 
numerical values of the quantities that characterize the MF behaviour. 
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2. The Model and the Multifractal Measure 

Long flexible polymer chains [|| provide a classical example of fractal objects (see 
fig. |T]a). Of particular interest is the manner in which such a chain absorbs diffusing 




Figure 1. a) A SAW of N steps each of size a in d-dimensional space. In the 
asymptotic limit N — > oo the mean square distance between the extremities is 
related to the 'mass' iV by A ~ (R/a) 1 ^ with v{d = 2) = 3/4, v{d = 3) = 0.588. 
This defines the fractal dimension of a SAW as D = 4/3, D = 1.70 in two and 
three space dimensions, b) A star of n RW (trajectories of diffusing particles) that 
reach the same point close to the absorbing polymer. The averaged nth moment 
of concentration of diffusing particles p scales as < p n >~ (R/a)~ Xn . The set 
of scaling dimensions A n depends on n in a non-trivial way, leading to the non- 
trivial behaviour of the flux moments of diffusing particles on the surface of fractal 
(polymer) absorber. 

particles. The flux of particles on the surface of the absorbing chain is defined by 
the density of incoming particles close to the absorber (see fig. |I|b). Taken p(r) to be 
the probability density of incoming particles the problem is to solve the steady-state 
diffusion equation (Laplace equation) 

V 2 p{r) = (1) 

with boundary conditions for a given absorbing polymer 

p(r) = on the surface of the absorber 

p(r) = poo = const for |r| — ► oo. (2) 

One motivation for the study of this problem is to gain insight to the description 
of diffusion limited aggregation (DLA) [ [TIJ. The last phenomenon is much more 
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complicated because of the fact that for DLA the boundary conditions are given 
on the surface of growing aggregate itself. The process described by equations ([[]), 
(0) may be rather considered as diffusion limited catalysis, when particles of one 
type interact with a prescribed fractal (a catalysing polymer) and transform into the 



other type [ O. 



Let the absorbing polymer of size R be chosen from the well defined ensemble of 
SAW. Then we are interested in the moments of the field p(r) close to the surface of 
the absorber (see fig. |l]b). For distances R » r > a (a being a cut off) the averaged 
moments < p{r) n > are expected to scale as [ |5j 

< p(r) n >~ (R/r)- x <- n \ (3) 

where < . . . > denotes the average over the ensemble of polymers. Let us take that 
the flux <f>(x) onto any randomly chosen point x of the absorber is proportional to 
the field p(a) at a point that is as close as a cut off length a from the absorber. 
Then, for the averaged moments of the flux one finds: 

< n >~ {R/a)- x(n) . (4) 

Following [ [5[] we associate with the flux of the Laplacian field a "harmonic measure" 
p{x) defined on the set of absorbing sites x : 

Ai(s)=0(aO/X>(*), ( 5 ) 

where the sum in (§) spans all sites x of the absorber. In particular, the mass M of 
the absorber is given by J2 X 1 = M and ([5]) may be rewritten as 

ti(x)=<t>(x)/{M <</>>), (6) 

where < . . . > denotes the site average < . . . >= XX- • \ 0- 

It is standard [ [3| to describe the properties of the harmonic measure by a set of 
exponents D(n) given by: 

5>(a;) n = (R/a)^ n)D{n \ (7) 

Note, that for n = D(0) is the fractal dimension of the absorber. Multifractal 
scaling is found when the spectrum D(n) is non-trivial: D(n) ^ -D(O), or corre- 
spondingly (c.f. eq. (H)) A(n) ^ nA(l). Note as well that from ([|), (^) the exponents 
A(n) and D(n) are related: 

nA(l) - A(n) = (l-n)[D(n) - D(0)]. (8) 

In many cases the description of a multifractal measure is performed in terms of the 
Holder exponent a and the spectral function /(a) that are defined by a Legendre 
transform of D(n) [ |3J]: 

a(n) = ±[( n -l)D(n)\, (9) 

f(a) = na(n) + {l-n)D(n). (10) 

The function f{a) differs for different fractal measures but it possesses certain 'uni- 
versal' properties. In particular, its maximum max a f{a) gives the fractal dimension 
of the absorber and / is convex /"(a) < 0. 
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3. The Multifractal Spectrum Exponents and the Spectral Func- 
tion 

The central idea that allows us to calculate the properties of the solutions of 
equation (ffl) with the boundary conditions (|2[) is that a path integral representation 
both for the field p(r) and for the absorber (being a RW or a SAW) is possible [ 
[|, |I2|| . In terms of the path integral solution of the Laplace equation one finds that 
p(r) at point r near absorber is proportional to the number of RW that end at point 
r and avoid the absorber. The nth power of this field is proportional to the nth 
power of the above mentioned number, i.e. it is defined by a partition function of 
a star polymer with n arms [ [13| (the latter is shown by dashed lines in fig. |l]b). 
Furthermore, introducing the mutual avoidance conditions between the "n-arm star" 
and the 2-arm polymer (representing the absorber) one has to calculate the partition 
function of a co-polymer star consisting of chains of two different species that avoid 
each other. These correspond to the trajectories of diffusing particles (being RW) 
and the absorbing polymer (which for the purpose of present study is chosen as RW 
or SAW). Making use of the previously developed [ |8| theory of copolymer stars and 
networks and mapping the model of co-polymer stars to the appropriate field theory 
[ [15| one may relate [ |T7| the MF spectrum exponents (|J) to the exponents that 



define the scaling properties of co-polymer stars. In particular, for a co-polymer 
star consisting of n\ chains of species 1 and n,2 chains of species 2 the number of 
configurations Z* scales like [||: 

~ (iz/a)^-™-^^ ( H ) 

Here, the rj nin . 2 represent the co-polymer star exponents. The latter have been calcu- 
lated using a field theoretical renormalization group approach [ ^|, [7j and are known 
in the third order of perturbation theory in e = 4 — d and pseudo-e [ [H| expansions 



[ D, [16[]. By means of a short-chain expansion [ [19| the set of exponents t] nin2 can be 



related to the exponents A(n) (^), that govern the scaling of the nth moment of the 
flux onto an absorbing linear chain. Considering the absorber to be either a RW or 
a SAW we define the exponents: 

\ R w{nf = -vL (12) 

^SAw(n) U = ~V2n + V20, (13) 

here, we keep notations of [ || [16| for the expressions for exponents rj defined in 
different fixed points of the renormalization group transformation. Based on these 
expressions we get: 

X (A ™ ^-1)^ , n(n-l)(-l+n + 3C(3)) g 3 

arw{£) — en ^ 1 g , l i4 J 

3en n(7n - 18)e 2 
asaw{£) - — H ^ 



n(-149 - 42n + 108n 2 - 276C(3) + 540C(3)n)e 3 

2048 ' 



(15) 
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Table 1. Exponents \Rw{ n ) and XsAw{n) obtained in e and in pseudo-e expansion 
(A(r)) techniques. 



n 




WW 




^SAW 


1 


0.99 


0.99 


0.71 


0.71 


2 


1.77 


1.81 


1.31 


1.33 


3 


2.45 


2.53 


1.86 


1.92 


4 


3.01 


3.17 


2.34 


2.44 


5 


3.51 


3.75 


2.78 


2.94 


6 


3.95 


4.28 


3.19 


3.41 



Here, e — 4 — d, d is the space dimension, and ((3) ~ 1.202 is the Riemann zeta 
function. 

Using the perturbation expansions for the A exponents and the relations for A(n) 
and the spectral function some algebra leads to the corresponding expansions for a n 
and f{a n ): 



a RW {e) 



2 + 



{-2n + l)e 2 



-An - 3C(3) + 3n 2 + 6((3)n + l)e 3 



(16) 



e 2 n 2 



+ 



(2n 3 + 3((3)n 2 - 2n 2 )e 2 
8 



;i7) 



0-SAW\£) 



e (7 
2 "4 + - 



36n)e 2 



128 



(-149 - 276C(3) - 84n + 324n 2 + 1080C(3)n)e 3 

2048 '' 



e (18n 2 + ll)6 2 
fsAw(e) - 2-~ — + 

(216n 3 + 540C(3)n 2 - 42n 2 - 83 + 264C(3))£ 3 
2048 



(19) 



The exponents A(n) as well as the Holder exponents a and the spectral functions 
f(a) in terms of the pseudo-e expansion [ [L8| are given elsewhere [ [17], pOfl . 

As is well known, the series of type (|I"4^-(|T9"|), as they occur in field theory appear 
to be of asymptotic nature with zero radius of convergence. However, knowledge of 
the asymptotic behavior of the series as derived from the renormalization group 
theory allows us to evaluate these asymptotic series (see e.g. [ 0). To this end sev- 
eral procedures are available depending on the additional information known for 
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(Fig.2a) 




(Fig.2b) a 

Figure 2. The spectral function f(a) for absorption on (a) RW and (b) SAW . 
Solid curves: 1 - [2/1] Pade approximant for e 3 results, 2- [2/1] Pade approximant 
for pseudo-e 3 results; dashed curves: 3 - e 2 results without resummation, 4 - pseudo- 
e 2 results without resummation; stars - resummed e 3 results; boxes - resummed 
pseudo-e 3 results. 
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the series to be resummed. We use the Borel summation technique improved by a 
conformal mapping procedure that has served as a powerful tool in field theoretic 
calculations (see [ ^1 f° r example). 

Table [I] contains the results for the exponents \rw{ u ) an d ^sAw( n ) obtained in 
e and in pseudo-e expansion techniques at d = 3 from the corresponding values of 
exponents rjf lt f 2 [ |3| with the application of this resummation procedure. 

4. Conclusions 

We have calculated the spectral function that describes the scaling of the mo- 
ments of measure defined by diffusion near an absorbing polymer. These moments 
were calculated as averages over all configurations of the absorber instead of perform- 
ing a site average. Thus the interpretation of f(a) itself does not directly correspond 
to the standard picture, where f{a) is interpreted as the fractal dimension of a par- 
ticular subset of the measure [|||]. Nevertheless we derive f{o) in the same way from 
the spectrum of scaling exponents of the moments of a measure. As we will see below 
in the general properties of f{ot) also hold for the present definition. 

Our numerical results for the spectral function for space dimension d = 3 are 
presented in Figs. ||]a,|2|b. They were obtained from the series for a n and f(a n ) as 
functions of n. We show the results of the resummation procedure described above 
applied to the series in e- and pseudo-e- expansion. For comparison we also show the 
curve for direct summation of the e and pseudo-e series to the 2nd order (here we 
recover the second-order results of [||). In addition we have performed an analytical 
continuation of our series in the form of a [2/1] Pade approximant for the third order 
series. It is obvious that the direct summation of the last series fails to converge and 
gives reasonable values for a n , f(a n ) only for small values of n, i.e. near the maximum 
of f{u) at n — 0. The symmetry of the Pade approximant holds only in the region 
shown and may be an artifact of the method. On the left wing, where it coincides 
with the resummed results the Pade approximant provides a continuation that is 
compatible with the estimation for the minimal a value ct m i n = d — 2. The Pade 
result is a m m(s) = 1.333, a m i n (r) = 1.017 for the RW absorber and a min (e) = 1.250, 
«min( r ) = 1-013 for the SAW absorber, which is calculated here only from 3rd order 
perturbation theory. 

As can be extrapolated from the Pade approximant and as was shown also on 
the basis approximations for higher moments n [ ||, f{o) as it is defined here, will 
become negative near a min and a max . For this reason the identification of f(a) as 
the fractal dimension of some identifiable subset is not possible here. Also the ex- 
trapolation of the resummed data seems to indicate such a behavior. Note, however, 
that the perturbative approach, even in combination with resummation and ana- 
lytical continuation is still not capable to give reliable results for high values of the 
expansion parameters. In particular this method is only good near the maximum of 

/(«) [HI- 

Though the results obtained for a n and f(a n ) for a specific value of n differ in 
both approaches, the same curve f(a) is described with better coincidence for the 
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left wing of the curves, corresponding to positive n. The resummation techniques 
we apply have proven to be a powerful tool already in the field theoretical approach 
to critical phenomena and have lead to high precision values for critical exponents. 
We believe that the application of these methods to calculations of MF phenomena 
also allows to improve the reliability and comparability of these results. 

A generalization of our present approach to the study of spectral functions that 
describe the absorption at the core of an absorbing polymer star with a given number 
of arms will be published elsewhere. 

It is our pleasure to acknowledge discussions with Lothar Schafer, Bertrand Du- 
plantier, Alexander Olemskoi. We are indebted to Gleb Oshanin for attracting our 
attention to the references [ 12 1. Part of this work was supported by DFG (SFB 237) 
and the Minerva foundation. 
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